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ABSTRACT 

We study a specific correlation in spherical harmonic multipolc domain for cosmic microwave back- 
ground (CMB) analysis. This group of correlation between A£ = An, n = 1, 2 ... is caused by symmetric 
signal in the Galactic coordinate system. An estimator targeting such correlation therefore helps remove 
the localized bright point-like sources in the Galactic plane and the strong diffused component down to 
the CMB level. We use 3 toy models to illustrate the significance of these correlations and apply this 
estimator on some derived CMB maps with foreground residuals. In addition, we show that our proposed 
estimator significantly damp the phase correlations caused by Galactic foregrounds. This investigation 
provides the understanding of mode correlations caused by Galactic foregrounds, which is useful for 
paving the way for foreground cleaning methods for the CMB. 



1. INTRODUCTION 

Separation of the cosmic microwave background (CMB) 
signal from extragalactic and Galactic foregrounds (GF) is 
one of the most challenging problems for all the CMB ex- 
periments, including the ongoing NASA WMAP and the 
upcoming ESA Planck mission. The GF produces the ma- 
jor (in amplitude) signal in the raw maps, which is local- 
ized at a rather small latitude band b < 30°. To avoid any 
contribution of the GF to the derived CMB map, start- 
ing from COBE to WMAP experiments, a set of masks 
and disjoint regions of the map are in use for extraction 
of the CMB anisotropy power spectrum (Bennett et al. 
2003a, b,c; Hinshaw et al. 2003b; Tegmark, de Oliveira- 
Costa & Hamilton 2003; Eriksen et al. 2004). The ques- 
tion is, what kind of assumption about the properties of 
the foregrounds should we apply for the data processing 
and what criteria determines the shape and area of the 
mask and the model of the foregrounds? To answer these 
questions we need to know the statistical properties of the 
GF to determine the strategy of the CMB signal extraction 
from the observational data sets. 

These questions are even more pressing for the CMB 
polarization. Unlike temperature anisotropies, our knowl- 
edge about the polarized foregrounds is still considerably 
poor. Additionally, we have yet to obtain a reasonable 
truly whole-sky CMB anisotropy maps for statistical anal- 
ysis, while obtaining a whole-sky polarization map seems 
to be a more ambitious task. Modeling the properties of 
the foregrounds thus needs to be done for achieving the 
main goals of the Planck mission: to the CMB anisotropy 
and polarization signals for the whole sky with unprece- 
dented angular resolution and sensitivity. 

Apart from modeling the foregrounds, Tegmark, de 
Oliveira-Costa & Hamilton (2003) (hereafter TOH) pro- 
pose the "blind" method for separation of the CMB 
anisotropy from the foreground signal. Their method 



(see also Tegmark & Efstathiou (1996)) is based on min- 
imizing the variance of the CMB plus foreground signal 
with multipole-dependent weighting coefficients w(£) on 
WMAP K to W bands, using 12 disjoint regions of the 
sky. It leads to their Foreground Cleaned Map (FCM), 
which seems to be clean from most foreground contam- 
ination, and the Wiener-Filtered Map (WFM), in which 
the instrumental noise is reduced by Wiener filtration. It 
also provides an opportunity to derive the maps for com- 
bined foregrounds (synchrotron, free-free and dust emis- 
sions . . . etc.). Both FCM and WFM show certain levels of 
non-Gaussianity (Chiang et al. 2003; Bershadskii & Sreeni- 
vasan 2004; Schwarz et al. 2004), which can be related to 
the residuals of the GF (Naselsky et al. 2003). Therefore, 
we believe that it is imperative to develop and refine the 
"blind" methods for the Planck mission, not only for bet- 
ter foreground separation in the anisotropy maps, but also 
to pave the way for separating CMB polarization from the 
foregrounds. 

The development of "blind" methods for foreground 
cleaning can be performed in two ways: one is to clarify 
the multipole and frequency dependency of various fore- 
ground components, including possible spinning dust, for 
high multipole range and at the Planck High Frequency 
Instrument (HFI) frequency range. The other requires 
additional information about morphology of the angu- 
lar distribution of the foregrounds, including the knowl- 
edge about their statistical properties in order to con- 
struct realistic high-resolution model of the observable 
Planck foregrounds. Since the morphology of the CMB 
and foregrounds is closely related to the phases (Chiang 
2001) of aj> >m coefficients from spherical harmonic expan- 
sion AT(f?, 4>), this problem can be re- formulated in terms 
of analysis of phases of the CMB and foregrounds, includ- 
ing their statistical properties (Chiang & Coles 2000; Chi- 
ang, Coles & Naselsky 2002; Chiang, Naselsky & Coles 
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2002; Coles ct al. 2004; Naselsky et al. 2003; Naselsky, 
Doroshkevich & Verkhodanov 2003, 2004). 

In Naselsky & Novikov (2005), it is reported that a ma- 
jor part of the GF produces a specific correlation in spheri- 
cal harmonic multipolc domain at A£ = 4: between modes 
ae^ m and cn + ^ m . The series of 4n-correlation from the GF 
requires more investigation. This paper is thus devoted to 
further analysis of the statistical properties of the phases 
of the WMAP foregrounds for such correlation. We con- 
centrate on the question as to what the reason is for the An 
correlation in the WMAP data, and can such correlation 
help us to determine the properties of the foregrounds, in 
order to separate them from the CMB anisotropies. 

In this paper we develop the idea proposed by Nasel- 
sky & Novikov (2005) and demonstrate the pronounced 
symmetry of the GF (in Galactic system of coordinates) 
is the main cause of the 4n correlation. The estimator 
designed in Naselsky & Novikov (2005) to illustrate and 
tackle such correlation can help us understand GF man- 
ifestation in the harmonic domain, leading to the devel- 
opment of "blind" method for foreground cleaning. In 
combination with multi-frequency technique proposed in 
Tegmark & Efstathiou (1996); Tegmark, de Oliveira-Costa 
& Hamilton (2003), the removal of An correlation of phases 
can be easily used as an effective method of determination 
of the CMB power spectrum without Galactic mask and 
disjoint regions for the WMAP data. It can serve as a 
complementary method to the Internal Linear Combina- 
tion method (Bennett et al. 2003c; Eriksen et al. 2004) 
and to the TOH method as well, in order to decrease the 
contamination of the GF in the derived maps. Such kind 
of correlation should be observed by the Planck mission 
and will help us to understand the properties of the GF in 
details, as it can play a role as an additional test for the 
foreground models for the Planck mission. 

This paper is organized as follows. In Section 2 we de- 
scribe the df m estimator for 4n-correlation in the coeffi- 
cients ae, m and its manifestation in the observed signals. 
In Section 3 we apply the estimator on 3 toy models which 
mimics Galactic foregrounds to investigate the cause of 
such correlation. In Section 4 we discuss the connection 
between the 4n correlation and the WMAP foreground 
symmetry. We also examine the power spectrum of the 
estimator and the correlations of df m estimator in Sec- 
tion 5. The conclusion is in Section 6. 



coordinate system, Yi^ m (9, 4>) are the spherical harmonics, 
\at, m \ and &e, m are the amplitudes (moduli) and phases of 
£, m harmonics. The superscript A in df' characterizes 
the shift of the £-mode in df m and, following Naselsky & 
Novikov (2005), we concentrate on the series of correlation 
for A = An, n = 1,2,3.... Note that the singal of Galaxy 
mostly lies close to 9 — 7r/2-plane. The estimator df m in 
form Eq.(l) is closely related with phases of the multipoles 
of the AT (9, (j)) signal. Taking Eq.(2) into account, we get 



dt, m — 



(3) 



From Eq.(3) one can see that, if the phase difference 
*£+A,m - ®e,m -» 0, then 

dtm - ae, m e^ /2 sin($ £+A , m - <f>t, m ) -» 0. (4) 

If $£+A,m — ®e,m <C 7r/2, the map synthesized from the 
dfm estimator is simply a map from the a^ iTO with phases 
rotated by an angle ±7r/2 and the amplitudes lessened 
by a factor | sin($£ + A, m — &t,m)\, while for non-correlated 
phases $£+A,m, &t,m we have specific (but known) modu- 
lation of the a^ m coefficients (see the Appendix). 

A non-trivial aspect of df m estimator is that it signifi- 
cantly decrease the brightest part of the Galaxy image in 
the WMAP K-W maps. In the following analysis we use 
a particular case n = 1 so that A = 4, although it can be 
demonstrated that for n = 2,3... the results of analysis 
do not change significantly as long as A < i?noiso, where 
£ no ise is the multipole number in the spectrum where the 
instrumental noise starts dominating over the GF signal. 



2.2. The An correlation in the WMAP data 

In this section we show how the df m estimator trans- 
forms the GF image in the WMAP K-W maps, taking 
from the NASA LAMBDA archive (?). In Fig.l we plot the 
maps synthesized from the df m estimator for WMAP K- 
W band (for A = 4 and £ max = 512) 



2. THE An CORRELATION AND ITS MANIFESTATION IN 
THE WMAP DATA 

2.1. The df m estimator 

As is shown in Naselsky & Novikov (2005) to illustrate 
the 4n-correlation, we recap the estimator taken from the 
combination of the spherical harmonic coefficients ae, m , 



at-, 



\ae,; 



(1) 



\at+A,m\ 

where |m| < £, and the coefficients a^ m = 
\at, m \ exp(z<I>^ m ) are defined by the standard way: 

AT{9,^) = Y j \at, m \exp(i*i,m)Ye, m (6,4>). (2) 

1=1 m=-l 

AT(9, 4>) is the whole-sky anisotropies at each frequency 
band, 9, (j) are the polar and azimuthal angles of the polar 



D(9^) = ^T J2 4,mYt, m (0,<t>). 



(5) 



1=1 m= 



Note that the amplitudes are significantly reduced in each 
map. It should be emphasized that the D(9, <j>) map is not 
temperature anisotropy map, as the phases are altered. 

Let us discuss some of the properties of the df m esti- 
mator, which determine the morphology of the D(9, 4>) 
maps. First of all, from Eq.(l) one can find that for 
all m = modes, the estimator is equivalent to zero if 
sign(a^o) = sign(a£ + A,o) an d it is non-zero (and doubled) 
if sign(a^o) = — sign(a£ + A,o)- In terms of phase difference 
in Eq.(4) this means that for m = modes df m estimator 
removes those which have the same phases, while doubles 
the amplitudes of others whose phases differing by 7r in 
the D{6, <p) maps. 




Fig. 1. — The D(9, 0) maps (from the top to the bottom) synthe- 
sized from the estimator for WMAP K, Ka,Q, V and W maps. 
Note that the colorbar limits from top to bottom are [—0.50,0.50], 
[-0.50,0.50], [-0.50,0.50], [-0.40,0.40] and [-0.38,0.50], respec- 
tively. 

However, such specific case of df , m estimator for m = 
modes is not unique for A — An only. It seems typical 
for any values of A parameter. What is unique in the 
WMAP data is that for A — An the order of sign for 
m = modes leads to the D(9, <fi) image without strong 
signal from the Galactic plane. 

We present in Fig. 2 the images synthesized of the even 
and odd m modes from the WMAP W band signal. The 
even and odd m modes reflect different symmetry of the 
signal, related to the properties of the spherical harmonics 
and the corresponding symmetries of the foregrounds. For 
even m the brightest part of the signal is mainly localized 
in the Galactic plane area (the top panel), while for odd 
m modes the signal has less dominated central part from 
the GF, but it has well presented periodic structure in 9 
direction (horizontal stripes) , from the north to south pole 
caps crossing the Galactic plane. In Fig. 3 we present the 
symmetry of the GF for the W band signal for even and 
odd I harmonics, including all corresponding m modes. 
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Fig. 2. — Comparison between even and odd m modes. The top 
panel is the map synthesized of the even m modes for 2 < I < 512 
from the WMAP W band signal. The 2nd is the D{9, <f>) map syn- 
thesized from the d^ m estimator on the 1st panel. The 3rd is that 
synthesized of the odd m modes from the WMAP W band signal, 
and the 4th is the D(8, </>) map on the 3rd panel. All the maps are 
plotted with colorbar limit [—0.5, 0.5] mK. 

As one can see from Fig. 2-3, the even m and the even 
£ maps (the top of Fig. 2 and Fig. 3) have a common sym- 
metrical central part, which looks like a thin belt covered 
in 9 ~ 7r/2 area and all < <f> < 2tt range. For odd modes 
the brightest GF mainly concentrate locally in 9 and <f> 
rectangular areas. Additionally, for the maps of even and 
odd t harmonics in Fig. 3 we have periodic structure of the 
signal in 9 direction, which is determined by the proper- 
ties of the spherical harmonics, and, more importantly, by 
the properties of a^, m coefficients of decomposition, which 
reflect directly corresponding symmetry of the GF. 
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Fig. 3. — Comparison between even and odd £ modes. The top 
panel is the map synthesized of the even £ modes for 2 < £ < 512 
from the WMAP W band signal. The 2nd is the D(8, <f>) map syn- 
thesized from the d^ m estimator on the 1st panel. The third is 
that synthesized of the odd £ modes from the W band and the 4th 
the D(9,<j>) map on the 3rd panel. All the maps are plotted with 
colorbar limit [—0.5, 0.5] mK. 

3. WHY DOES THE 4N CORRELATION APPEAR? 

In this section we want to examine why the An correla- 
tion appears. In order to answer this question we intro- 
duce 3 toy models for the Galaxy emissivity, which reflect 
directly different symmetries of the Galactic signal. In 
Appendix we analyze more general situation. These 3 toy 
models are the belt, the rectangular and the spots models. 
All 3 models are the simple geometrical shapes added with 
the WMAP ILC map. 

1. The belt model: 

We add on top of the WMAP ILC map with 

AT b (0,4>) = A = const, (6) 

if tt/2 - S < 9 < tt/2 + <5,and 4> C [0, 2tt], where S is 
the halfwidth of the belt. We set A — 10 mK and 
6 = 5°. 

2. The rectangular model: 

We add on top of the WMAP ILC map with 

AT r (9,<j>) =A = const (7) 

for tt/2- 6 < 9 < tt/2 + 5, and 7 < 4> < 2tt-j. We 
further use two variants 



• 6 = 5° and 7 = tt, 

• (5 = 5° and 7 = 7r/4. 

3. The spots model: 

This model is to mimic the properties of symmetric 
bright point-like sources convolved with a Gaussian 
beam B(9, 4>) with FWHM=5°. So on ILC map we 
add 



AT S (9,0) = J2 / A ] 6 D (cos9 3 )5 D ( ( f>-<p ] ) 
j J 

xB(9-9 j} <p- fa^jdfa , (8) 

where 5 D is the Dirac-5 function, and the ampli- 
tudes of the point sources Aj are in order of 10 
mK. 

Below we examine these toy models to see how they can 
illustrate the An correlation. 

3.1. The belt model 

Ideologically this model best illustrates that the symme- 
try of the GF signal can help remove the GF itself without 
any additional assumption about properties of the fore- 
grounds and of Galactic mask as well. The theoretical 
basis is that the properties of the phases for the belt sig- 
nal are related to m — modes (see below). Thus, any 
method that can remove all m — modes from the map is 
automatically able to remove the belt-like Galactic signal 
and reproduce the ILC signal for all m ^ modes with- 
out any errors. However, these m = modes contribute 
significantly to the low multipole part of the power spec- 
trum Ce for the reconstructed CMB signal (which is in 
our model the ILC signal for £, m 7^ modes). Using df m 
estimator we can avoid the problem of reconstruction of 
the m = modes and the corresponding power spectrum 
of the D(9, <j>) map (see Section 5). 

In Fig. 4 we plot the resultant D(9, 4>) maps for the belt- 
like model. As one can see, the particular case seems pe- 
culiar for the d^ m estimator, which does not remove the 
Galactic signal properly. According to the private com- 
munication with Eriksen et al., this simple model of the 
Galactic signal should reflect the properties of spherical 
harmonics, namely the An correlation mentioned in Nasel- 
sky & Novikov (2005). The symmetry of this model is 
related to 2n correlation of even harmonics, while for odd 
harmonics the Galactic signal must vanish. However, for 
some of the m = harmonics the df m=0 estimator differs 
from zero, and the corresponding stripes determines the 
morphology of the image (see Fig. 4, the second from the 
top). To show that effect clearly, in Fig. 4 (the third) we 
plot the map, which contains the non-zero df m=(j modes. 
Remarkably, the morphology of defects in the second and 
third one are practically the same. To demonstrate that, 
in the D{9, <fi) map (the second in Fig. 4) we simply remove 
all non-zero dg m _ q modes and get the map shown in Fig. 4 
as the 4th. So, the second and the 4th maps are the same, 
excluding all non-zero modes df m=0 from the D(9, 0) map. 
This result is not surprising, if one takes into account the 
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symmetry of the belt signal, all the phases oil, m harmon- 
ics for m ^ are exactly those of t, m 7^ harmonics in the 
ILC map, while all a^ m= o in the first from the top in Fig. 4 
map are mainly determined by the belt signal (Galaxy). 
To show that in Fig. 4 we simply remove all ai^ m= o modes 
from the ILC plus the belt map and the resultant is shown 
in Fig. 4 the bottom map. One can see, that this is the ILC 
map for all ae_ rn ^o harmonics. What properties of this toy 
model are crucial in producing the An multipolc correla- 
tion of the phases? Taking the definition of the belt signal 
into account, one can obtain that for the belt signal 




2A [5 e , e=0 - P e 1 (COS 6) Sin @Si t 2n$m,m=o] , (9) 




where n = 1,2,..., S x , y is the Kronecker symbol, 
F^ 1 (cos 9) is the Legendre polynomials and 9 = tt/2 — 5. 
For > 1 the asymptotic of a\ m is (Gradshteyn & 
Ryzhik 2000) 



AAT(£) 



e,m=o 



2kT(£+§ 



sin 



'2n- 



(10) 



Then for dg t o from Eq.(2) we obtain 



d A - 
a e,o — 



AAT(£) 



| sin [(£ + ±)6] | 
\sin[(l + A + ±)5] 







|sin 





■ sm 



(l + A + -)S 



(11) 



As is seen from Eq.(ll), dg depends on the sign of the 
first and second terms in the brackets. Let us assume that 
for some I — I* the first term is positive and the second 
term is negative. In this case 



Fig. 4.— The belt model. The 1st panel is the WMAP ILC 
map plus artificial uniform Galactic belt signal (see text). The 2nd 
panel is the D(9, 0) map synthesized from d^ m -estimator (£ < 250) 
on the 1st panel. The 3rd panel is also the D(8, <f>) map from djt - 
estimator (the same as the 2nd), but with all m = modes in 
being set zero. The bottom panel is synthesized from the ai m of 
the 1st panel, but with all m = modes in the ai m being set zero. 
The colorbar limit for the first panel is [—0.6, 10.5]mK, and the rest 
arc [-0.5,0.5]mK. 

Thus, for a given "belt" model of the Galaxy emissivity 
for the range of multipoles £5 ^> 1 the 4n multipole corre- 
lation does not exist. However, in the opposite approxima- 
tion £5 <C 1 the properties of the multipole coefficients are 
determined by the asymptotic of Legendre polynomials in 
Eq.(9) 



sm 



sm 



A + -)S 



(12) 



and unlike compensation of the modes df Q = 0, we have 

the signal in order of magnitude close to — 4A/y/2~7r£~i . 
Moreover, taking into account that the amplitude of that 
signal decreases very rapidly in comparison with ampli- 
tudes of the CMB signal, starting from some multipoles 
£ = £ cr r~j S^ 1 the order of sign in Eq.(2) is determined 
by the CMB signal and the number of horizontal stripes 
increase rapidly, as it is seen from the Fig. 4 (the second 
map from the top). 



4,; 



T(£ + 2) 



2 I 



r(¥ 



(13) 



T(£)T (1 

and formally a\ m — for all £ = 2n, n = 1,2..., but 
£6 <C 1. As follows from Eq.(13), the symmetry of the 
belt model does not require 4n correlation at all. If A pa- 
rameter is even A = 2n, the sum £ + A is even too, and 
the corresponding a b lrn coefficient vanishes. The case for 
A = An is a special one for the more general correlation 
A = 2n, which will be broken at the limit £8 ^> 1. How- 
ever, the belt model as discussed above provides a useful 
estimation of the properties of the d^ m estimator for more 
general cases, when the symmetry of the GF is not so high 
as in the belt toy model. Particularly, to prevent any con- 
tribution to the D(6, 4>) map from the highly symmetrical 
part of the GF, we will use further generalization of the 
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df' m estimator as follows: 



a e,n 



: at,; 
0, 



if m = 0. (14) 



We will use this generalization for df m in all the subse- 
quent discussions (excluding descriptions of the toy models 
in §3.2 and 3.3). 

3.2. The rectangular model 

Let us discuss the properties of the model which mimics 
the Galaxy image in the map as a rectangular area, char- 
acterized by halfwidth 8 in 9 direction and $ in <fi direction 
and with a constant amplitude of the signal A in the area 
vr/2 - 5 < 9 < vr/2 + S, $ < cj) < 2vr - $. For that model 
the corresponding ae_ rn coefficients are 



-2A 



(2£ + 1)T(£ - m + 1) 



4-KT(£ + m + 1) 

*7r/2+r5 

: / dO sin 9P^ (cos 9). 

'tt/2-6 



sin(m<i>) 



(15) 



The properties of the integral in Eq.(15) depend on the 
parameter IS. If IS < 1, then Pp(cos9) ~ P £ m (0), where 



pp(0) = J^T 



r i 



r 



1 — £ — to 



while for M > 1 we get (Gradshteyn & Ryzhik 2000) 

2 r(^ + m + l) 



(16) 



Pf(cos^) 



\/27r sin( 



r(£ + |) 



2^* 



X cos 

Thus, for these asymptotics we have 
T(£-m+l) 



3 - 

2 ' 

1 7T 

— run 

2 4 



r(*- 

if £<5 < 1, and 

T(. 



i)r(i + ^2)r(i 



(17) 



(18) 



1) 



T(£ + m+ 1) 

r .tt^ + v 

x ^ COS( — 



1) 



r(* + g) 



(2£+l)g 
tt(£ + m) 



sm 



7r(£ + m) 



1 



(19) 

if » 1. For odd £ + to, as seen from Eq.(18), a r t m = 

if £ S -C l 3 . That means that the main term in Eq.(15), 
proportional to 5 /it is related with the even harmonics 
£ + m = 2n, n = 1, 2 The leading term, which deter- 
mines the sign of a£ +A m is r( 1 ~ £ ~ 2 "'~ A ) j n the denomi- 
nator of Eq.(19). For that term we have 



r( 



l-l-m- A, 



' l-l-m^ 



(20) 



■J 



and for A = 4k, k = 1, 2 . . . the sign of a£ +A rn is the same, 
as for a r lrn . Taking into account that for d\ m estimator 
the order of signs for £, m and £ + A, m is crucial, we can 
conclude that the compensation of the central part of the 
signal requires A = 4k, k = 1,2.... However, it does not 
guarantee that for £5 ^> I modes A = 4fc criteria leads the 
compensation of the signal. To show this let us describe 
the asymptotic £8 ^> 1, when the symmetry of the a^ m 
coefficients is determined by Eq.(17). As one can see from 
this equation, if £ + m = 4k, k = 1,2... the sign of a[ +A m 

is now determined by the combination sin [(£+ h)S] and 
does not show any 4fc-correlation at all. Moreover, ac- 
cording to the properties of the sine mode the shift of the 
argument £ by the factor A = 4k just transforms it to the 
combination 

i ' i 

cos(4kS) 
sin(4kS). (21) 





Fig. 5. — The rectangular model. The 1st panel is the WMAP ILC 
map plus the rectangular signal (see text). The 2nd is the D(8,4>) 
map synthesized from the d^ m -estimator without Eq. (14) general- 
ization (£ < 250) on the 1st panel. The 3rd is from the 1st panel, 
but with all m = modes in the a? m being set zero. The bottom 
is the D(d ! (f>) map from the d^ m -estimator (without Eq.(14) gen- 
eralization) on the ILC map alone. The colorbar limit for the 1st 
panel is [—0.6, 10.5]mK, and the rest are [— 0.5, 0.5]mK. 

Thus, one can see that 4fc-correlation requires some re- 
striction on the 5-parameter 

4kS = 27rm. (22) 
Thus, for k = 1 and to — 1 the halfwidth of the rect- 
angular area must be close to S = tt/2. If, for example, 



The first non- vanished for odd £ + m is in order of magnitude {IS) 
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S = ir/A then we will have correlation for k — 2, but not for 
k = 1. Practically speaking, for ^ > 1 we can have some 
particular symmetry, but not general symmetry A = An. 
Conclusions concerning this rectangular model of GF are 
clearly seen in Fig. 5. 

3.3. The spot model 

To understand how each sort of defects is related with 
corresponding 4n-correlation, we introduce the model of 
defects, which can be describe as a sum of peaks with 
amplitudes Aj and coordinates 8j , <f>j . For analytical de- 
scription of the model we neglect the beam convolution of 
the image of point sources (PS), but we include it in the 
numerical simulation. For the model of defects 

AT(0, <t>) = ^ A^icosO - cos 6 j)8{(j) - fa). (23) 
i 

For the a^ m coefficients of the spherical harmonics expan- 
sion from Eq.(6) we get 

As for the rectangular model, we will assume that all 
9j = ir/2 and simply we will have ai^ m ~ -P™(0) for 
18 1, as well as 18 <C 1. As it was shown in pre- 
vious section, P™(0) clearly demonstrate An correlation. 
We would like to point out that for the spots model this 
correlation now is strong, unlike the belt and rectangu- 
lar models. 4 Moreover, implementation of the Gaussian 
shape of the PS which come from beam convolution does 
not change that symmetry at all. To show that, in Fig. 6 
we plot the model of two PS with amplitudes in order to 
10 mK, combined with the ILC map. The reason for such 
effect is quite obvious. The beam convolution does not 
change the symmetry of the model, but rescale the am- 
plitudes of the PS by a factor exp[— £(£ + l)]/2c 2 , if we 
assume the Gaussian shape of the beam. 




Fig. 6. — The D(9,(f>) maps for PS model. The top panel is the 
ILC map plus 2 symmetrical PS signal and the bottom is the D(8, 0) 
map (£ < 50) without Eq. (14) generalization. The colorbar limits 
are [— 0.6, 6.3]mK and [— 0.5, 0.5]mK, respectively. 

An important question is that is the symmetry of the 
Galaxy image in </> direction important for extraction 
of the brightest part of the signal, or is the effect sim- 
ply determined by the symmetry of the Galaxy image 
in 9 direction? To answer this question, in Fig. 7 we 

4 Strongly speaking, we should use the terms I + m = 4n-correlation. 



plot the result from the df' estimation in the model 
with 5 spots located at 9j — tt/2 with different ampli- 
tudes and different <f)j. As one can see no symmetry in 
4> direction was assumed. The result of reconstruction 
clearly shows that location of the sources in the Galactic 
plane in (f> direction is crucial. Unlike the model with 
symmetric location of the spots in Fig. 6, now the resid- 
uals of the extraction of the spots dominate over the 
rest of the signal in the Galactic plane. However, as 
is seen from Fig. 7, the 4n-correlation of image exists. 




Fig. 7. — The D(9, 0) maps for another PS model. The top is the 
ILC map plus 6 PS signal and the bottom is the D(8, 0) maps from 
the top panel without Eq. (14) generalization. The colorbar limits 
are [— 0.6, 6.3]mK and [— 0.5, 0.5]mK, respectively 

The properties of the ae. m coefficients in the spots model 
are related with the sum (see Eq.(24)) 

5(m) = Y^Aje-™^. (25) 
i 

Actually, Eq.(25) determines the phases of the a,£ im coeffi- 
cients. Let us discuss the model of two symmetrically situ- 
ated PS with the same amplitudes A± — A 2 and fa = tt/2, 
fa = 2tv — <fii. For this particular case Im[iS(m)] = 0, while 
Re[S(m)} = 2Acos(irm/2) . 

For m = 2k + l,fc = 0,1,2... we get Re[5(m)] = 
lm[S(m)] = and the contribution of the strong signal 
to the map ILC + two PS vanishes. This means that for 
amplitudes a^ m — ci >m +pe, m , where ct >m and p£ iTn corre- 
spond to the ILC and PS signals respectively) and phases 
^e.rn of ai^ rn coefficients, we have 

0<e,m=2k+l = Ce,m=2k+1, 

^,m = 6,m=2fe+l, fc = 0,l,2... (26) 

where are the ILC phases. As one can see, this is 
a particular example, when strong, but symmetric in <f> 
direction signal do not contribute to the set of a^ >m coeffi- 
cients at least for the defined range of multipoles. 

Let us discuss the other opposite model, in which the 
number of spots in the galactic plane is no fewer than 
2, and their <f>j coordinates are random in some range 
$, 2tt — No specific assumptions about the amplitudes 
are needed. In this model the sum S(m) in Eq.(25) mostly 
is represented by m = modes, 5(0) = J2j Aj, while all 
m/0 modes S(m) -C 5(0) because of randomness of the 
phases mfa. This model, actually is close to the rectan- 
gular model, in which the width of rectangular side in the 
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4> direction now is <&, 2tt — <&. At the end of this section 
we would like to demonstrate, how the symmetry of the 
Galaxy image in <f> direction can determine the properties 
of the D(8, (j> map. For that we rotate the W-band map by 
20° along the pole axis and produce the same estimation 
of df m , as is done for the Galactic reference system. The 
result of estimation is shown in Fig. 8. For comparison, in 
this Figure we plot the difference and sum between D(6, <f>) 
maps before and after rotation. From these Figures, these 
new symmetry of the W band map after rotation simply 
increase the amplitude of signal in Galactic plane zone, 
especially in the the central part of it. 

At the end of this section, we summarize the main re- 
sults of investigation of the given models of the GF signal. 




For highly symmetrical signal, like the belt model, 
all ai, m coefficients vanish for the multipole num- 
bers I = 2n, n = 1, 2 . . ., but 15 < 1. Modification 
of the df m estimator in form of Eq.(14) is crucial to 
prevent any contribution to the D(9, </>) map from 
the GF. 



•,£«• 7: '■<•■ if'" '-'. 



Less symmetric model, like the rectangular model, 
requires An transition for the multipole numbers for 
df m estimator, which appears for the range of mul- 
tipoles ft < 1. If the resolution of the map we 
are dealing with is low, max 5 < 1 that An correla- 
tion appears for all £ < £ max and the corresponding 
a>t,m coefficients for GF are in order of magnitude 
ai,m ~ A(£5) 3 . For 15 » 1 the An correlation of 
phases does not exist at all. 




Fig. 8. — The W band map, rotated by 20° in <j> direction (top). 
The 2nd is the D(8, 0) map derived from the top panel (without 
Eq.(14) generalization). The 3rd is the difference between D(9,4>) 
maps from W band before and after rotation. The bottom is the 
sum between D(8, cf>) maps (before and after rotation). The colorbar 
limits are [— 0.5,0. 5]mK for all panels. 



• The amplitude of the GF signal, A, and its de- 
pendency on I (like A{1) = A)(^/£ m in)~' 3 , 4nin is 
the minimal multipole number for which A(£ m i n ) 
achieve the maxima, and (3 is the power index) 
are crucial for establishing of the An correlation 
of phases. Taking asymptotic 15 <C 1 into ac- 
count, and defining the critical multipole number 
£ cr r~j <5 -1 , we can estimate the corresponding am- 
plitudes A(£ cr ) = y4(l/<5£ m in) - ' 3 . If at that range of 
£ we get A(£ cr ) < C*{1), where C{£) is the CMB 
power spectrum, the An correlation would be estab- 
lished for all range of multipoles I < £ maX) even if 
it vanishes for the GF signal for 15 ^> 1. Start- 
ing from t = £ cr and for £ > £ cr the corresponding 
ai tin for GF play a small role in amplitude noise, in 
comparison to the amplitudes of the CMB signal. 



• The df m estimator effectively decreases the am- 
plitudes of the point-like sources located in the 
Galactic plane, if they have nearly the same ampli- 
tudes and are symmetrically distributed in <fi direc- 
tion around Galactic center. Non-symmetrical and 
different in amplitudes point-like sources after im- 
plementation of df m estimator produce significant 
residues. 



4. SYMMETRY OF THE WMAP FOREGROUNDS 

In this section we apply the proposed df m estimator 
to the maps for WMAP Q, V and W band foregrounds 
(which are sum of synchrotron, free-free and dust emis- 
sion). We then transform them by the df m estimator. 
These foreground maps do not contain the CMB signal and 
instrumental noise, therefore they allow us to estimate the 
properties of the GF in details. In Fig. 9 we plot the D(8, 4>) 
maps for Q, V and W band foregrounds (A = 4) for the 
multipole range I < 46. This range is determined by the 
resolution of the WMAP foregrounds maps (£ < 50). As 
one can see from these maps, the GF perfectly follows to 
An multipole correlation, which remove the brightest part 
of the signal down to the level ± 50 mK for the Q band, 
-0.19,0.50 mK for the V band , -0.09,0.29 mK for the 
W band and —0.1, 0.1 mK for the D(9, cf>) map, the differ- 
ence between V and W foregrounds. Note that these limits 
are related with the brightest positive and negative spots 
(point sources) in the maps, while diffuse components have 
significantly smaller amplitudes. To show the high resolu- 
tion D{6, (p) map which characterizes the properties of the 
foregrounds in V and W band, in Fig. 10 we plot the map 
of difference V — W bands, and the corresponding D(9, 4>) 
map for I < 250. Note that V — W map does not contain 
the CMB signal, but for high I the properties of the signal 
are determined by the instrumental noise. 
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5. THE POWER SPECTRUM AND CORRELATIONS OF THE 

d MAP 

5.1. What is constructed from the df' m estimator? 

To characterize the power spectrum of the D(8, <p) maps 
we introduce the definition 





D(£) = 



-A 



2£ + 1 ^ 



(27) 



If the derived df m signal is Gaussian, that power repre- 
sents all the statistical properties of the signal. For non- 
Gaussian signal, D(£) power characterizes the diagonal el- 
ements of the correlation matrix. From Fig. 11 it can be 
clearly seen that for WMAP foregrounds, especially for V 
and W bands, the power spectra of D(8, (f>) are signifi- 
cantly smaller than the power of the CMB, for estimation 
of which we simply use the power of TOH FCM map, 
transformed by df m estimator as 




Fig. 9. — (From top to bottom the D(8, (f>) map for Q, V and 
W band foregrounds, respectively, with f max = 46, A = 4. The 
bottom is the D(8, </>) map from the V and W band map difference: 
V -W. 



D fcm (£) = 



21- 



1 ^ 

m— — j 



R,. 



\Ce+A,m\ 



| 2 , (28) 



assuming that FCM map is fairly clean from the fore- 
ground signal. An important point of analysis of the 
WMAP foregrounds is that for V and W bands df m esti- 
mator decreases significantly the amplitude of GF, practi- 
cally by 1 to 2 order of magnitude below the CMB level. 

The most intriguing question related to 4n-correlation 
of the derived map from the WMAP V and W band sig- 
nals is what is reproduced by the df m estimator? The 
next question, which we would like to discuss is why the 
power spectrum of df m estimation of the V and W bands 
shown in Fig. 12 are practically the same at the range of 
multipoles I < 100, when we can neglect the contribution 
from instrumental noise to both channels and differences 
of the antenna beams. The equivalence of the powers for 
these two signals, shown in Fig. 12, clearly tell us that these 
derived maps are related with pure CMB signal (which we 
assume to be frequency independent). 




Fig. 10. — The 1st panel is the difference between V and W bands: 
V— W, the 2nd is the same map, but with all m = modes in m 
being set zero, and the 3rd is the D(6, </>) map from the 1st panel. 
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10 

multipole (1) 



Fig. 11— The power spectrum G{£) = D(£)£(£+l) /2w for D(6, 0) 
map (W band, the solid line with stars) in comparison with the G(£) 
power of the D(8,<f>) map for FCM map (the thick solid line). The 
dash line represents the power spectrum G(£) for the D(8, (p) map 
for the difference between Q and W bands: Q— W, and the dash-dot 
line is for V— W. 

In this section we present some analytical calculations 
which clearly demonstrate what kind of combinations be- 
tween amplitudes and phases of the CMB signal in the 
V, W bands and phases of foregrounds are represented in 
the df m estimator. As was mentioned in Section 1, this 
estimator is designed as a linear estimator of the phase 
difference &e+A,m — &e,m, if the phase difference is small. 
Let us introduce the model of the signal at each band 

a< itm = c £,™ + ^im> where- Q ;TO is frequency independent 
CMB signal and Ff^ is the sum over all kinds of fore- 
grounds for each band j (synchrotron, free-free, dust emis- 
sion etc.). 




10 100 
multipole (1) 




ur" I , , , , , , 

10 

multipole (1) 

Fig. 12.— The power spectrum G(£) = D(£)£(£ + l)/2ir for D 
map (W band, stars) in comparison to the powers of the ILC map 
(dash) and the FCM (solid line). The dotted line represents the G 
power of the D map for difference V— W. The bottom is the low 
resolution power spectra D{£) for W and V bands. 

According to the investigation above on the foreground 
models, it is realized that without the ILC signal the df m 
estimation of the foregrounds, especially for V and W 
bands, corresponds to the signal 5 

dfi f) = F l>m - \pA_ Fi+A , m , (29) 

the power of which is significantly smaller then that of the 

5 Hereafter we omit the mark of channel j to simplify the formulas 



CMB 

jAAcmb) \C£m\ /on\ 

d e „ = C Lm - |C£ +A ,m- (30) 

\Ce+A,m\ 

In terms of moduli and phases of the foregrounds at each 
frequency band 

Fe, m = \Fe,m\exp{i®e,m), 

ce, m = \ct m \ exp(i& im ), (31) 

where &e,m and are the phases of foreground and the 
CMB, respectively. And from Eq.(31) we get 

4£ f) = \F t>m \ (e***-™ - e**<+*"» ) , (32) 

and practically speaking, we have $£ >m = <&e+A.rn ■ Thus, 
taking the 4n correlation into account, we can conclude 
that it reflects directly the high correlation of the phases 
of the foregrounds, determined by the GF. Moreover, if 
any foreground cleaned CMB maps derived from different 
methods display the 4n correlation of phases, it would be 
evident that foreground residuals still determine the sta- 
tistical properties of the derived signal. 



5.2. An phase correlation of the D map 

One of the basic ideas for comparison of phases of two 
signals is to define the following trigonometric moments 
for the phases m and $?t,m as: 

1 £ 

Cs(£,e') = ^ Yl cos (&', ro -*/, m ); 

v *• m=l 
1 1 

v 1 m=l 

(33) 

where I < i . We apply these trigonometric moments 
to investigate the phase correlations for TOH FCM and 
WFM. For that we simply substitute I = i' in Eq.(33), 
and define ^ >TO as the phase of FCM and ^e,m as that of 
WFM. The result of the calculations is presented in Fig. 13. 

From Fig. 13 it can be clearly seen that the FCM has 
strong A£ — 4 correlations starting from £ ~ 40 which 
rapidly increase for I > 40, while for WFM these correla- 
tions are significantly damped, especially at low multipole 
range I < 40. However, the df m estimator allow us to clar- 
ify the properties of phase correlations for low multipole 
range. The idea is to apply df m estimator to FCM and 
WFM, and to compare the power spectra of the signals 
obtained before and after that. According to the defini- 
tion of df m estimator, the power spectrum of the signal is 
given by Eq.(28), which now has the form 

D W = jYj l Q > m l 2 t 1 - cos (&+A,™ - tt, m )\ ■ (34) 

m 
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To evaluate qualitatively the range of possible non- 
Gaussianity of the FCM and WFM, in Fig. 14 we plot the 
function F{£) = 2[D(£) - 2C{tj\/[D{£) + C{tj\ for FCM 
and WFM, in which we mark the limits ±2/^/1. As one 
can see, potentially dangerous range of low multipoles is 
£ = 3,4, £ = 21 -24, £ ~ 100 - 150 for the WFM. Non- 
randomness on some of the multipole modes is mentioned 
in Chiang & Naselsky (2004). 
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Fig. 13. — The Cs(£) and Si(^) trigonometric moments for the 
cross-correlation of phases between the TOH FCM and WFM (the 
first pair). The solid line represents the limit when the phases are 
identical. The middle pair is The Cs(£) and Si(£) trigonometric 
moments for the FCM with phase difference t;e+4,,m — ^t,m- The 
bottom pair is Cs(^) and Si(^) trigonometric moments for the WFM 
with phase difference ^+4 m ~ £,1 m- 

The last term in Eq.(34) corresponds to the cross- 
correlation between £, m and £ + 4,m modes, which should 
vanish for Gaussian random signals after averaging over 
the realization. For a single realization of the random 
Gaussian process this term is non-zero because of the same 
reason, as well known "cosmic variance" , implemented for 
estimation of the errors of the power spectrum estimation 
(see Naselsky et al. 2004). Thus 



D{£) 



-V 

£ ^ 



and error of D(£) is in order to 

AD(£) 2 



D{£) 



(35) 



(36) 



Fig. 14.— F(t) function for the TOH FCM (top) and the WFM 
(bottom). The dotted lines represent ±2/ \f~t limit. 

At the end of this section we would like to demon- 
strate that application of df~ m estimator to maps with 
foregrounds residuals, such as the FCM, provides addi- 
tional "cleaning". In Fig. 15 we present the Cs{£) and 
Si(^) trigonometric moments for the FCM with shift of 
the multipoles £ = £ + 2. One can see that the A = 2 
correlation of phases is strong (practically, they are at the 
same level as A = 4 correlations). However, after df m 
filtration these correlations are significantly decreased. 

The implementation of the df^ estimator to the non- 
Gaussian signal significantly decreases these correlations. 

The properties of the df m estimator described can mani- 
fest themselves more clearly in terms of images of the CMB 
signal. In Fig. 16 we plot the results of the maps with df m 
implemented on FCM and WFM, in order to demonstrate 
how the estimator works on the non-Gaussian tails of the 
derived CMB maps. In Fig. 16 we can clearly see that the 
morphology of the D(8, (f>) maps are the same and differ- 
ence between Df cm (9, <j>) and D w f m (6, </>) is related to point 
sources residuals localizes outside the galactic plane (see 
the 3rd panel). A direct substraction of the WFM from the 
FCM reveals significant contamination of the GF residuals 
and non -galactic point sources ( the third from the bottom 
and bottom maps) . The second from the bottom map cor- 
responds to difference between Df cm (9,4>) and .D w f m (0, 0) 
for which the amplitudes of the signal represented in col- 
orbar limit ±0.1 mK. One can see that the GF is removed 
down to the noise level. In combination of the phase anal- 
ysis we can conclude that the implementation of the d^ m 
estimator looks promising as an additional cleaning of the 
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GF residuals and can help investigate the statistical prop- 
erties of derived CMB signals in more detailed. 
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Fig. 15. — The Cs(£) and Si(£) trigonometric moments for the 
FCM at A£ = 2 (the first pair). The 2nd pair is for At = 2 after 
df= 4 filtration. The 3rd and 4th pairs are for Al = 1 before and 
after df , ~ i filtration, respectively. 



Fig. 16. — The comparison between the FCM and WFM be- 
fore and after implementation of the d^ m estimator. From the top 

to the bottom: -Df cm (6> ,<f>) map for the FCM, _D w f m (9,0) for the 
WFM, the difference between Df cm (0, <j>) and D w f m (0, </>), and the 
difference between FCM and WFM. All these maps are plotted with 
colorbar limit ±0.5 mK. The last pair are the same as previous pair, 
but for the colorbar limit ±0.1 mK. For all the ITlcLpS {-max — 500. 



6. CONCLUSION 

In this paper we examine a specific group of correla- 
tions between £, which is used as an estimation of the sta- 
tistical properties of the foregrounds in the WMAP maps. 
These correlations, in particular, among phases are closely 
related to symmetry of the GF (in Galactic coordinate 
system). An important point of analysis is that for the 
foregrounds the correlations of phases for the total fore- 
grounds at V and W bands have specific shape when 
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^e,m — &i+A.m, A = An, n = 1,2,3.... These correlations 
can be clearly seen in the W band of the WMAP data 
sets down to £ max = 512 and must be taken into ac- 
count for modeling of the foreground properties for the 
upcoming Planck mission. We apply the df m estimator to 
the TOH FCM, which contains strong residuals from the 
GF and show that these residuals are removed from the 
D w { m (6,4>) map. Moreover, in that map the statistics of 
the phases display the Gaussian statistics closer than the 
original FCM (no correlation of phases between different 
l,m modes except between I + A, to and I, m, which is 
chosen as a basic one, defined by the form of df m estima- 
tor.) 

In this paper we do not describe in details the properties 
of the signal derived by df m estimator from the WMAP V 
and W bands. Further developments of the method, 
including multi-frequency combination of the maps and 



CMB extraction by the estimator will be in a separate pa- 
per. To avoid misunderstanding and confusion, here we 
stress again that any D(9, fa maps synthesized from the 
df m are by no means the CMB signals (since the phases of 
the these signals are not the phases of true CMB) and the 
true CMB can be obtained after multi-frequency analysis, 
which is the subject of our forthcomming paper. 
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APPENDIX 



In this Appendix we would like to describe general properties of the An periodicity of the Galactic signal, taking into 
account its symmetry. We adopt the following model of the signal, which seems to be general. Let define some area 
around Galactic plane S — X^jLi s pix ~ Ns p i X , where s pix is the pixel area and index j mark the location of the pixel. We 
assume for simplicity that all the pixels in the map are have the same area. In polar system of coordinates corresponding 
angles 6j and fa mark the position of j-th pixel in the map. Let us define the amplitude of the signal per each pixel as 
Tj . Thus the map which corresponds to the Galactic signal is now 

AT(6,fa = ^T.,c5(cos0 - cos6j)6(<t> - fa). (Al) 

3 

Let assumes that Galaxy image is localized in ^-direction as 7r/2 — 5 < 6j < tt/2 + 5 and it could be or could not be 
localized in 0-direction. Additionally we will assume that signal per each pixel Tj is the sum of Galactic foreground signal 

Tj and CMB plus instrumental noise signal Tj. Important to note that statistical properties of these two components 
are different as in terms of amplitudes, as in terms of pixel-pixel correlations (TjTk). Particularly, in the area S we have 
Tj ^> Tj, while outside S we assume that Tj <C Tj. Using proposed model of the signal in the map we can obtain 
corresponding at. m coefficients of the spherical harmonic expansion 



a t , m = ^|^| E W(cos^e-^ , (A2) 

which can be represented as a sum of foreground F^_ m coefficients and the CMB plus noise coefficients C£_ rn . In order 
to understand the nature of 4n-periodicity of the Galactic foreground, let discuss the model when q ;TO = O.Then, from 
Eq.(l) the subject of interest would be the phases of foreground $£ m related to the Fe, m coefficients as follows 

J2 i Tj P™ (cos 6j ) sm(mfa ) 

tan$i iTO = - J — — — , (A3) 

l^j Tj Pp(cos 0j) cos(m<pj) 

where sum over j corresponds to the pixel in the area S. Let's define the difference of phases, using their tangents. 

N 

tan $ e +A,m - tan $ e ,m = , (A4) 

where 

= E TjTk(P™ (cos 0j)Pil A (cos 9k) sinm^ cosm<pk — P7 1 ( cos ^fe)-f > i+A( cos ^i) sinm^fc cosm^; 
D = J2 T,T k P^ l (cos 9,)P^ A (cos 9 k ) cos mfa cos mfa. 

(A5) 

As one can see from Eq.(A4), if 3>^ TO ~ &£+A,m, then tan$^ jm — tan$^ + A, TO — ®e,m — ®e+A,m , which determine 
the properties of df m estimator for correlated phases (see Eq.(l). Below the object of out investigation is function N 
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from Eq.(A5). Particularly we are interesting in asymptotic N — > 0, which should reflect directly the symmetry of the 
foreground signal. Simple algebra allows us to represent N function in the following form 

N = Y J TjT k {PF{cos6 3 )PF +A {cose k ) - P^cos 9 k )P% A (cos 9 3 )) sinm(fo - cj>j). (A6) 

Taking into account that area S is located close to the 9 = ir/2, let us discuss the properties of N-function at the limit 
£tt/2^> 1, using asymptotic of the Legendre polynomials. After simple algebra we obtain 

N -zH-r^T^W rtPM 3WTA Z A)smm(<fa - <j>j), A7 

2G% A (9j,6 k ) = {cos[(^ + + fc ) + tott - |] + cos[(^ + ^)(6j - 9 k )]}(cos9 k A - cos^A) 

-{sin[(* + ±)(0j + 6 k ) + m-K - |] - sin[(f + - ^)]}(sin# fe A - sin^A). (A8) 

From Eq.(A8) one can see the symmetry of the Legendre polynomials which manifest themselfs trough cos 6^ A and svnOjA 
modes. If 6 3 - = 7r/2, then depending on A we will have 

cos(^A) = 1, A = 4n,n = 1,2...; 



where 



2 

.TT 

2' 



sin(^A) = 0, A = 2n,n = 1,2.... 



(A9) 

Thus, choosing A = 4n, n = 1,2... mode we take the corresponding properties of the Legendre polynomials into consid- 
eration. However, as one can see from Eq.(A9) A = 4n,n = 1,2... periodicity of the Galactic image is not exact. In 
reality we have pixels which containt galaxy signal having thetcij coordinate close to tt/2, but not exactly equivalent to 
it/2. Let us introduce a new variable Sj — tt/2 — 9j,5j <C 1 characterized the deviation of the j-th pixel location from 
the theta = tt/2 plane. From Eq.(A8) one can find that for A = 4n,n = 1,2... cos^A ~ cosinSj ~ 1 — (4n^)2/2. 
Thus , if pixel j containt the signal from the galactic foreground, the deviation from the center of the galactic plane 
should be small enough:^ <C 7r/4n. It is clear that this condition does not necessarily correspond to the properties of 
the Galactic image, which is clearly seen from the K, Ka and Q band signals. Taking the above-mentioned properties 
of G™ A=4n (^ : j, 9 k ) function, we represent the asymptotic of this function at the limit 5j <C 7r/4n, which is applicable for 
analysis of the Galactic signal at V and W bands. 

2G£ A=4n (0,-A) = A2 %~ — [cos[(e+ l -){5 k 5,)] + (-l)^cos^ + \){5 k + 5,)] 

+A(S 3 S k ) |sin[(^ + \){5 k - 5 3 )\ + (-lf +m sm\{t + \){5 k + *,-)]} . (A10) 

Thus, combining Eq.(A7) and Eq.(AlO), we obtain 

N TxY,TjTkGT A =in{6j,6 k )smm{<f> k - <j> 3 ). (All) 

j,k 

One may think that choose of described above 4n-mode of Legendre polynomials automatically guarantee cancellation 
of the brightest part of the signal from the map without any restriction on symmetry and amplitude of the foreground. 
To show that symmetry of the Galactic signal is important, let us discuss a few particular cases, which illuminate this 
problem more clearly. 

Firstly, let take a look at galactic center (GC), which is one of the brightest sources of the signal. For the GC 
corresponding amplitudes T 3 are localized per pixels, for which 4> 3 ~ in the galactic system of coordinate. From 
Eq.(All) one can see that for GC the function N is equivalent to zero . More accurately, taking into account that image 
of the GC has characteristic sizes 86 = 5<j) ~ FWHM, where FWHM is the Full Width of Half Maximum of the beam, in 
Eq.(All) additionally to 6 3 <C 1 parameter we get small parameter mFWHM -c 1. 

Secondly, let discuss the model of two bright point like sources, located symmetrically relatively to the GC. Let assumes 
that for that point sources T\ ^ T 2 , but 2 — 4>i — tt. Once again, from Eq.(All) we get N = for all m and these point 
sources will be automatically removed by the df^ estimator even if they have <5i 7^ 62 ■ 

Another possibility related to the symmetry of the Galactic image in 9 direction. We would like to remind, that 
Eq.(All) was obtained under approximation £6j <C 1, where 0j = tt/2 — Sj. This means, that 15 j can be as big enough 
(£5j > 1), as small (15 j < 1) as well. For 15 j < 1 from Eq.(AlO) we obtain 

1G? A = in (ej,6 k ) = A2( ^~ Sl \ l + (-l) e+m ] + A(£+ ±)[-(5j S k f + (-lf +m (5 2 j Si)]. (A12) 

As one can see from Eq.(A12) the bright sources located on the same 9 coordinates (Sj = 5 k ) docs not contribute to dj n m 
estimator. 
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